Potential for ill-posedness in several 2nd-order formulations of the Einstein equations 
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Second-order formulations of the 3+1 Einstein equations obtained by eliminating the extrinsic 
curvature in terms of the time derivative of the metric are examined with the aim of establishing 
whether they are well posed, in cases of somewhat wide interest, such as ADM, BSSN and generalized 
Einstein-Christoffel. The criterion for well-posedness of second-order systems employed is due to 
Kreiss and Ortiz. By this criterion, none of the three cases are strongly hyperbolic, but some of them 
are weakly hyperbolic, which means that they may yet be well posed but only under very restrictive 
conditions for the terms of order lower than second in the equations (which are not studied here). 
As a result, intuitive transferences of the property of well-posedness from first-order reductions of 
the Einstein equations to their originating second-order versions are unwarranted if not false. 

PACS numbers: 04.20.Ex, 04.25. Dm 



I. INTRODUCTION 

A common practice to study the well-posedness of a 
time-dependent second-order system of partial differen- 
tial equations is to reduce the system to first order. This 
is done because first-order systems of PDE's have been 
amply examined and are currently very well understood. 
By reducing to first-order, what is meant is to find a first- 
order system of PDE's whose solution space contains a 
subset that is equivalent to the solutions of the origi- 
nal second-order system by a trivial identification. If the 
first-order reduction is well posed, then the solutions of 
the first-order system are bounded in terms of the initial 
data of the first-order system, and theorems of existence, 
uniqueness and stability under small perturbations then 
follow. 

It is tempting to presume that the well-posedness of 
a first-order reduction is a sign of well-posedness in the 
originating second-order system, especially because the 
originating second-order system is much smaller (in terms 
of the number of variables) and may be more amenable 
to a numerical implementation. However, in general rel- 
ativity, the existence of constraints leads to first-order 
reductions that are equivalent to the original second or- 
der problem with regards to constrained solutions only. 
The evolution equations of two first-order reductions of 
the Einstein equations that differ by linear combinations 
with the hamiltonian or momentum constraints are not 
equivalent. Thus, whether a well-posed first-order re- 
duction of the Einstein equations guarantees the well- 
posedness of the corresponding second-order problem is 
not at all clear. 

The quintessential model of a well-posed time- 
dependent second-order equation in three dimensions is 
the wave equation in flat space. By the well-posedness of 
the wave equation, more complicated systems of PDE's 



can be shown to be well posed if they consist of series of 
wave equations, as is the case of the 3+1 Einstein equa- 
tions in the harmonic gauge or time-harmonic gauge 
But beyond the wave equation, little is known about gen- 
eral well-posedness criteria for systems of second-order 
time dependent PDE's. 

A criterion for the well-posedness of time-dependent 
second-order systems of PDE's has been developed re- 
cently by Kreiss and Ortiz |2| ■ Even though the criterion 
applies to cases that are much more general than what 
we are interested in, for our purposes it can be stated 
as follows. Consider a system of n second-order partial 
differential equations with constant coefficients. The sys- 
tem is of the form 



jk 



(1) 
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where u is a vector containing the n fundamental vari- 
ables and, for each value of (j,k), Ajk is an n x n 
constant matrix. An overdot denotes a partial deriva- 
tive with respect to the time coordinate (d/dt), so that 
/ = d 2 f/dt 2 , and ,j= djdxK For any arbitrary unit 
covector £j = we define the n-dimensional matrix 

P Q (£) = J^^k A jk£j(,k- The system is strongly hyperbolic 
if and only if the eigenvalues of Po(£) & r e strictly posi- 
tive and Po(£) has a complete set of eigenvectors which 
is uniformly (in £) linearly independent. 

This criterion is entirely similar to strong hyperbolic- 
ity of first-order systems Q except for the condition that 
the eigenvalues must be strictly positive, which excludes 
the value zero as well as all negative values. The reason 
for positivity of the eigenvalues is, essentially, that each 
eigenvalue must allow for two waves travelling with the 
same speed but in opposite directions (so, in a sense, the 
eigenvalues are the squares of real characteristic speeds of 
any sign). The reason for excluding the vanishing eigen- 
value is that there is only one associated "wave" with 
zero speed instead of two. This criterion emulates fully 
the case of the three-dimensional wave equation. Any 
second-order system of PDE's that satisfies this criterion 
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is well posed in the standard sense, that is: its solutions 
are bounded by the initial data irrespective of the spec- 
tral frequency of the data 2\. On the other hand, if 
the eigenvalues are real and non- negative and/or there 
isn't a complete set of eigenvectors, the system is re- 
ferred to as weakly hyperbolic. According to @, weakly 
hyperbolic systems can develop "catastrophic exponen- 
tial growth" when adding lower order terms or consid- 
ering variable coefficients. It is inferred that if the lin- 
earization of a system of nonlinear second-order equa- 
tions around a constant background is weakly hyperbolic, 
then the associated non-linear system itself is prone to 
"catastrophic exponential growth" , as is its linearization 
around any background that is not constant. Examples of 
what Kreiss and Ortiz refer to as "catastrophic growth" 
appear in [2|. 

As far as we are aware of, the criterion has been devel- 
oped only for linear systems with constant coefficients, 
but there is reason to presume that it can be general- 
ized to variable coefficients and quasilinear systems in a 
manner similar to the case of first-order systems. There- 
fore, in the following, we will use this criterion to analyze 
the well-posedness of second-order formulations of the 
Einstein equations insofar as they are linearized around 
flat space. The case of the standard ADM equations 3 
is studied in Section [H] The case of the widely used 
BSSN equations is dealt with in Section InTl Finally, 
the case of the generalizeded Einstein-Christoffel equa- 
tions (EC) @, which includes the case of the Einstein- 
Christoffel equations themselves 0, is developed in Sec- 
tion IIVI We find that none of the second-order versions 
of the equations is strongly hyperbolic. The relevance of 
this result is summarized in Section Ivl 

II. THE ADM EQUATIONS IN SECOND 
ORDER FORM 

Throughout the article we assume the following form 
for the metric of spacetime g ao in coordinates x a — (x % , t) 
in terms of the three-metric 7^ of the slices at fixed value 
of t: 

ds 2 = -a 2 dt 2 + ~,,d.v'(l.i J , (2) 

where a is the lapse function. The Einstein equations 
G a b = for the four-dimensional metric are equivalently 
expressed in the ADM form [4|: 

7y = -2aKij, (3a) 
Kij = a (Rij - 2K a K l i + KKij) - D^a, (3b) 

with the constraints 

C = i (fl - + K 2 ) = 0, (4a) 

C, = DjK^i — DiK = 0, (4b) 

to be imposed on the initial data. Indices are raised with 
the inverse metric r ) lJ , Di is the covariant three-derivative 



consistent with 7^ , i£y is the Ricci curvature tensor of 
7y, R its Ricci scalar, Kij is the extrinsic curvature of 
the slice at fixed value of t and K = 7^ Kij . This system 
of evolution equations is a (partial) first-order reduction 
of the original second-order Einstein equations, which we 
can recover by substituting in (|3b|) the extrinsic curvature 
in terms of the time derivative of the metric as given by 

% = -2a 2 R lJ +j lm j a %u~l:l lmA fhn'it 3 + -it 3 +'2aD i Dja. 

2 a 

(5) 

Linearizing around flat space, so that 7y = Sij + hij and 
a = 1 + e, these equations read 

hij = S kl (hkiaj — hu^kj — hji^i + hij t ki) + 2e,ij . (6) 

Obviously the problem depends on the specification of 
the lapse function. We are here interested in two special 
cases: the case of constant unit lapse and the case of 
lapse equal to the square root of the determinant of the 
three-metric. 

In the first place, then, suppose e = 0. The equations 
then read 

hij = 5 kl {hki,ij — hu t kj — hjiki + hij^ki) (7) 

and conform to Eq. Q). Solving the eigenvalue problem 
of the principal symbol -Po(£) is equivalent to making the 
ansatz hij = Vij exp i(^a; fc — st) for an arbitrary covector 
£i, which yields the following: 

s 2 V l3 = 6 kL (V kl ^ 3 - Vu&Zj - Vjfoti + Vij&b) (8) 

or, equivalently 

s %. = + t&v - t,jV ik e - iiV jk e , (9) 

where V denotes the trace of Vij and indices are raised 
with <5* J . The eigenvalues are thus labeled by s 2 . One 
can easily see that this problem admits solutions with 
vanishing eigenvalue, as follows. Assume s = 0, which 
yields 

= + tejV - £jV ik Z k - ZiVj k £ k (10) 

Contracting with (5 IJ we find V — £ l Vy£ J = 0. Using this 
information back into (|10fl we have 

= + ti^v^r - ^v ik e - ZiV jk e (11) 

Now contracting with £ J yields 

= Zi(V-Z l VimD, (12) 

which are three identities. Therefore, three of the six 
components of Vij are free, the remaining three being 
given by (|llfl . One can pick the three free components 
to be the three projections Vij£? , in which case, by i|ll[l . 
the other three are vanishing. This observation leads 
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directly to the fact that there are three linearly indepen- 
dent eigenvectors with vanishing eigenvalue, the eigen- 
vectors being the six-dimensional unit vectors along the 
directions of the three projections Vij£ 3 . This is so for 
every arbitrary direction £j. As an illustration, rep- 
resenting the six-dimensional eigenvectors in the form 
u = {V xx ,V X y,V xz ,V yy ,V yz ,V zz ), for & = 6f = (1,0,0) 
we have that Vik£ = V xx , V xy , V xz are free, whereas, by 
(fTTfl . the remaining components V yy , V yz and V zz van- 
ish. An arbitrary eigenvector in the degenerate space of 
eigenvalue has thus the form 

°u = (V XX ,V XV ,V XZ , 0,0,0) 

= V XX {1, 0, 0, 0, 0, 0) + V xy (0, 1, 0, 0, 0, 0) 
+V XZ (0, 0,1, 0,0,0), 

where V xx , V xy and V xz are arbitrary real numbers. There 
are thus three zero-speed eigenvectors associated with the 
x— direction which can be chosen as 

V = (1,0,0,0,0,0), 
°u 2 = (0,1,0,0,0,0), 
°u 3 = (0,0,1,0,0,0), 

which are manifestly linearly independent. In accordance 
with the criterion of Section [IJ this is enough to demon- 
strate that the evolution equations @ are not strongly 
hyperbolic. Nonetheless, for completeness, one can cal- 
culate the remainder of the eigenvalues and eigenvectors 
in a similar manner (or one could, of course, proceed by 
any standard methods of linear algebra to calculate all 
the eigenvectors and eigenvalues at the same time). The 
reader can verify that the only other eigenvalue is s 2 = 1 , 
for which Eq. © reads 

o = e^v - ~ i,r,,r . (13) 

These six equations are all satisfied if and only if 

V^ k =^- (14) 

This in turn implies that the three components of Vij 
other than the projections V, k £, k are arbitrary. Denoting 
them by V-j they are given by 

v± = Vij - ^v ]k e - i,v lk i k + teifvu? (15) 

and are such that V^^ = by construction (which 
means that there are only three independent components 
in Vjj). This leads directly to the conclusion that there 
are three linearly independent eigenvectors with light 
speed, labelled by the three components of Vy other than 
Vik£ k , and they are, clearly, linearly independent from 
the zero-speed eigenvectors (labelled by Vi k ^ k ). For in- 
stance, in the case that = (1, 0, 0) as above, by ((151) the 
three arbitrary components of Vij are V yy ,V yz and V zz , 
which, if combined with l|14f> . leads to the fact that an ar- 
bitrary eigenvector in the three-dimensional degenerate 



space of eigenvalue 1 is of the form 

l U = (Vyy+V ZZ ,0,0,Vyy,Vy Z ,V ZZ ) 

= Vy y (l, 0, 0, 1, 0, 0) + V vz (0, 0, 0, 0, 1, 0) 
+^(1,0,0,0,0,1). 

where V vv ,V yz and V zz are completely arbitrary real 
numbers. This shows that three unit eigenvectors can 
be chosen as 

l Ul = -^=(1,0, 0,1, 0,0), 

x u 2 = (0,0,0,0,1,0), 
= -^=(1,0,0,0,0,1), 

which are manifestly linearly independent of each other 
and of °Ui. 

In summary, the principal symbol Po(£) of Q) admits a 
complete set of eigenvectors and has real but not strictly 
positive eigenvalues. The linearized second-order ADM 
equations with constant unit lapse are, thus, weakly hy- 
perbolic. Consequently, the original nonlinear second- 
order ADM equations JSJ are prone to catastrophic expo- 
nential growth, as is their linearization around any back- 
ground but flat. 

Suppose now that the lapse function is equal to the 
square root of the determinant of the three-metric. In 
the linearization, this means that 2e,ij — 5 kl h k i^j. With 
this choice of lapse, Eq. 10 reads 

hij — 5 kl (2hkiaj — hu^kj — hji^i + hij^i) (16) 
The associated eigenvalue problem is 

- vn + 2&jV - ijV lk e - ^v jk e (17) 

One can see that this eigenvalue problem admits two lin- 
early independent solutions for s = 0. To see this quickly 
(without necessarily using standard algebraic methods to 
solve the problem completely), set s — and contract 
with <5 ij , which yields V = (2/3)£%-£?. Substituting 
this back and contracting this time with £ 3 one has 

= l^ l VU m (18) 

Thus four out of the six equations are solved by setting 
CVij£ J = and consequently also V = 0. This implies 
that two components of Vij are free, which can be taken 
as the two Vyf? other than £*Vij£ J . The remaining two 
equations in the set fix the remaining two components of 
Vij in terms of these: 

= Vij-ZjVi k Z k -ZiV sk Z k (19) 

As in any eigenvalue problem, the fact that two compo- 
nents of Vij are left arbitrary by (|17|l with s = leads 
directly to the conclusion that there are two linearly inde- 
pendent eigenvectors (which can be calculated in any way 
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the reader finds appealing) associated with zero speed 
and any direction For instance, if £j = (1,0,0), then 
V xy and V xz are free, but V xx = by JTHJ, and V yv ,V vz 
and 14z also vanish by virtue of (|19|l . The generic eigen- 
vector in the degenerate space of s = is thus 

°u = (0,V xv ,V xz , 0,0,0) 

= 7 xy (0, 1, 0, 0, 0, 0) + V xz (0, 0, 1, 0, 0, 0), 

where V xy and V xz are arbitrary real numbers, so there 
are only two linearly independent eigenvectors and they 
can be chosen as 

V = (0,1,0,0,0,0), 
°u 2 = (0,0,1,0,0,0). 

This is enough to conclude that the linearized second- 
order ADM equations with "harmonic" lapse, namely 
Eq. I|ltj|) . are not strongly hyperbolic. 

For completeness, one can calculate the remaining 
eigenvalues and eigenvectors. The reader can verify that 
the only eigenvalue other than zero is s 2 = 1, which has 
three eigenvectors associated with it, namely the compo- 
nents of Vij other than Vik£ k - In keeping up with our 
illustration but skipping over the procedure, which is en- 
tirely similar to the one used three times in the proceed- 
ing, for £i = (1,0,0) the three eigenvectors associated 
with light speed are (i.e., can be chosen as) 

V = -^(1,1,1,0,0,0), 

l u 2 = (0,0,0,0,1,0), 
l u a = -^=(0,0,0,1,0,-1), 

which are manifestly linearly independent of each other 
and of °wi,°M2- Because the eigenvalues are real and 
non-negative and there is not a complete set of eigenvec- 
tors, Eqs. p6fl are weakly hyperbolic. One can infer that 
the second-order nonlinear ADM equations J5J) are prone 
to "catastrophic exponential growth" even in the case of 
a lapse function proportional to the determinant of the 
three-metric, as is their linearization around any nonflat 
background. 

This result does not contradict Q- In |2|, the authors 
do consider the linearized second-order ADM equations 
and conclude that in the case of a lapse function that is 
proportional to the determinant of the three metric the 
problem is strongly hyperbolic. However, they impose, 
additionally, the linearized constraints on the evolution 
equations, which, as a consequence, do not have the form 
(|16fl . Thus the problem that Kreiss and Ortiz found to be 
well posed is equivalent to a constrained evolution prob- 
lem for the ADM equations in second-order with "har- 
monic" lapse, where the solutions are constrained at ev- 
ery time slice. In contrast, we have demonstrated here 
that the corresponding unconstrained evolution problem 
is not well posed. 



III. THE BSSN EQUATIONS IN SECOND 
ORDER FORM 



In the case of vanishing shift vector, the formulation 
of the Einstein equations referred to as BSSN 5] consists 
of the following evolution equations 



= -2aA, 



6' 



-K 



K = DiDja + a(A t jA lj + -K 2 



(20a) 
(20b) 

(20c) 



2a 



-2A 1 - 



(20d) 



1 



+2T k \iT l)kl + T kl iT klj - 2DiD j( f, + 4Di(f>Dj 

-pi j (t k , k +f kli (2f ikl + f kH )) 

+%(&D l 4 - 2D 1 ^) - 
i a 

+KA lJ - 2A a A\ 



TF 



(20e) 



for the 15 variables 






= Y^ ln ( d et7*j) 


(21a) 


7ij 




(21b) 


K 


= l lJ K l0 


(21c) 


Aij 


= ^ {k* ~ paKj 


(21d) 


f' 


= "7*4 


(21e) 



As in the ADM case, the initial data for these evolution 
equations must be chosen to satisfy the constraints 
Other than as applied to the initial data, we will disre- 
gard the constraints. 

We obtain a six-dimensional second order system for 
the metric variables by substituting Ay in terms of 7y- 
back into (|20e|) . and K in terms of cb into (|20cl) . Ad- 
ditionally, we substitute T l back in terms of 7 y ,j into 
(j3Ue|) . 

The linearization around flat space 7y = Sij + h^ 
implies that 7y = Sij + and <fi — h/12, where 
h.^ = hij — (l/2)Sijh and h — 8 lJ hij. After linearization 
around flat space, the result is the following second-order 
system 

1 



<j> = -8 ij e Hj 



(22a) 



hij = 5 ( hij^i — hii^j — hji^ki + -SijS rs hi S} k, 



+4 



(22b) 
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Consider first the case of constant lapse equal to 1. 
The equations reduce to 



= 



8 [ hij^ki — hiLkj — hji t ki + -^5ij5 rs hi s ^r 



+4 



(<f>,ij - g fly <* W 



and the corresponding eigenvalue problem, with hij = 

Vij expi(^kX k — st) and <j> — -0expi(^x fc — st), is 



than (*Va£ l . There are two other eigenvalues different 
from zero. One is s 2 = 1, with three associated eigenvec- 
tors which have free values of ip and of the components 
(23a) of V i:j other than = 8^ip. The other one is negative 
s 2 = — 1 , and has one eigenvector associated with it which 
has : rj) = and % = (3/2)(&& - (l/3)<%)£ fe V fe ^ with 
£ Vjy£ arbitrary. So 1251) has a complete set of eigen- 
vectors and its eigenvalues are real but not non-negative, 
and thus it is not even weakly hyperbolic. The conse- 
quences to the full nonlinear BSSN equations in second- 
order form and with "harmonic" lapse are the same as in 
the case of unit lapse. 



(23b) 



s 2 ^ = 



(24a) 



s 2 % = % - SMI? - ZjVut 1 + -6ijeV k £ 



(24b) 



This system has a vanishing eigenvalue s 2 = with 
three linearly independent associated eigenvectors, la- 
belled by the three components Vyf . In order to see 
this, set s 2 = in l|24bj) and contract with £ l . This 
yields £j(i/j - (l/8)£ fe t4;£') = 0. So three equations 
are satisfied by the choice ip = (l/8)£ fc Vfcz£ z . Addition- 
ally, since s 2 = 0, then also (|24all is identically satis- 
fied. Four out of the six equations are thus satisfied with 
this one choice, which means that three of the fields are 
free. This is enough to conclude that the second-order 
BSSN equations with unit lapse, namely Eq. I|23|l . are 
not strongly hyperbolic. We can calculate the remain- 
ing eigenvalues and eigenvectors. The reader can verify 
that there are two other eigenvalues. One is s 2 = 1, 
with two associated eigenvectors, which are the compo- 
nents of Vj other than V«£ , with = and ip = 0. 
The other eigenvalue is negative s 2 = —1/3, and has 
one eigenvector associated with it, which is given by 
ip = 0, Vu£} — except for £ l Vu£ l which is free, and 
Vj = (3/2)(&& - (l/3)^)e fc Vfci?'. Thus there is a com- 
plete set of eigenvectors and all eigenvalues are real, but 
are not non-negative, and the equations are not even 
weakly hyperbolic. This means that the full nonlinear 
second-order BSSN equations with unit lapse have the 
potential for "catastrophic growth", as does their lin- 
earization around any background other than flat. 

Things change very little if one considers now the case 
of a = «/ det jij , or e,y = 6(p,ij. Equations l|22(l reduce 
to 



<P = S^cb,. 



(25a) 



5~ [ hij^ki — hu^kj — hji^ki + -5ij5 rs hi s ^kr 



+ \§(<j) n j-^8 l j8 kl <j) 1 ki)- 



(25b) 



This system has two eigenvectors with s 2 = 0, which are 

given by ?Va? = 0, ^ = and % = 

with arbirary values for the two components Vu^ other 



IV. THE GENERALIZED EC EQUATIONS IN 
SECOND ORDER FORM 



The generalized Einstein-Christoffel (EC) formulation 
is first introduced in [|| , where the method of derivation 
from the ADM equations is described without the details 
of the resulting equations themselves. The principal part 
of the evolution equations of the system appear explicitly 
in Hi as follows 



Hi 
Ktj = 

fkij — 



= -2aKi 



-a^OJk.j 
-adkKij + . 



(26a) 
(26b) 
(26c) 



where fkij constitute a set of 18 first-order variables de- 
fined by the following relationship with the first deriva- 
tives of the three-metric 

7y,fc = 2fkij+Vlk{i (fj)s s - f s j)s)+ 1 4 lij(fks s -f s ks)- 

(27) 

Here ij is a free parameter. The generalized EC family 
requires the lapse function to be densitized, that is: to be 
proportional to the square root of the determinant of the 
three-metric. The original (standard) EC formulation 
obtained by Anderson and York 7j corresponds to the 
choice of rj = 4. 

We are interested in the second-order system of equa- 
tions for 7jj that is implied by this 30-dimensional first- 
order problem. We start by inverting l|27|l in order to 
have an expression for fkij that we can use to substitute 
in terms of 7^ . The inversion yields 



2 fkij — 



7ij,k + lik{l U lls,j ~ l lS ll 3 ,s) 

+7jk(7 ls lis, t -l ls lu,s) 

+ ^-lij(l lS lls,k - l lS llk,s)- 



(28) 



If one uses (|2"%1) in its left-hand side, Eq. H26cfl reduces 
to a linear combination of the components of the momen- 
tum constraint Cj. Therefore, this equation is redundant 
for the problem that we are interested in. However, elim- 
inating fkij from the right-hand side of l|26bjl by means 
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of (H yields: 

Kij = a (Rij - 2KuK l j + KK l3 )-D l D ] a+ T ^-^a ll3 C. 

Zrj 

(29) 

One can see that, at rj = 4, the evolution equation for 
the extrinsic curvature in the standard EC formulation 
is an exact transcription of the ADM evolution equation 
without mixing of the constraints, whereas for any other 
value of the parameter 77, the Hamiltonian constraint is 
involved. This means that the second-order version of the 
standard EC is indistinguishable from the second-order 
version of the ADM equations, dealt with in Section ITT1 
Additionally, one can also show that for rj — 12/7, the 
second-order version of the generalized EC formulation 
is indistinguishable from the second-order version of the 
BSSN formulation, which is dealt with in Section II I II 
Yet, for 77 7^ 4, 12/7, the second-order version of the gen- 
eralized EC is a genuinely different problem which re- 
quires separate study. 

The second-order problem that we seek is found by 
substituting K ij in terms of jij in the left-hand side of 
(|29|l and writing explicitly the right-hand side in terms 
of Jij,kl- 

a~ 2 7ij = l kl HjM - l kl lu,kj - l U lk .a + ^l M lki,i 3 

+ ^l t3 (l kl l mS l m s,ki ~ l kl l ms l ml ,ks) 

+ ... (30) 

In the linearization around flat space (7^ — <5y + hij) 
we have 

hij = S^hijj-i — 5 kl hu,kj — 5 kl hkj,u + 25 kl hki,ij 

+^S ij (S M 5 ms h ms , kl - 5 kl 8 ms h mlM ),{Zl) 

which has the following eigenvalue problem in terms of 
= Vij expi(£ fc x fc - st): 

s% = V ij -V ik t%-V jk { k t i +2Vt i £j 

+ ^&ij{V-t*Vuti l ). (32) 

This problem has the following eigenvalues. First we 
have s 2 = 1 with three linearly independent eigenvectors 
which are the three components of Vij other than the 
projections Vy£ J (= ^V). Then we have s 2 = 2(77 — 2)/r) 
with one eigenvector given by Vy = 2r)V(2£i£j + <%(?7 — 
i)/2rj)/(7rj — 12) with arbitrary V, which is linearly in- 
dependent of the other three except at 77 = 4. Fi- 
nally we have s 2 = with two eigenvectors given by 

= tiV jk e + tjV lk e with v = evki? = 0. 

So there is a complete set of eigenvectors for all values 
of 77 7^ 4, 12/7. However, there are two eigenvectors with 
s 2 = for all 77 ^ 2,4,12/7. Furthermore, if 77 = 2 
then s 2 = has multiplicity 3, but if 77 < 2 then there 
is a negative eigenvalue. Collectively, this means that 
for no values of 77 is the second-order problem implied 



by the generalized EC formulation strongly hyperbolic. 
The second-order problem is weakly hyperbolic for 77 > 2, 
which includes the ADM case, and this is consistent with 
Section [D] Yet the second-order problem is not even 
weakly hyperbolic in the range 77 < 2 considered in 0, 
which includes the BSSN case, and this is consistent with 
Section ITTT1 



V. CONCLUDING REMARKS 

Summarizing, we have found that the second-order 
equations for the three-metric that arc implied by a 
large number of first-order reductions of the 3+1 Ein- 
stein equations whose evolution equations differ by the 
addition of different multiples of the Hamiltonian con- 
straint are potentially ill-posed irrespective of whether 
the lapse is constant or densitized. We say potentially 
only because the linearized argument that we use does in 
on way rule out the possibility that the lower-order terms 
that are present in the nonlinear equations may yet (for- 
tituously) result in a well-posed nonlinear problem. This 
should raise a flag for numerical efforts seeking formu- 
lations for general relativity with the smallest possible 
number of variables. 

On the other hand, studying the second-order equa- 
tions implied by formulations with such disparity of prop- 
erties as ADM, BSSN and generalized EC gives us yet an- 
other perspective for what it is about such formulations 
that characterizes their properties. 

In the case of the ADM equations and their associ- 
ated second-order version, we have verified a fact that 
has been suspected for a long time, starting with 0, 
where the second-order version itself is found unfit for 
standard hyperbolicity studies. Subsequently, many au- 
thors have concluded that the ADM equations themselves 
are ill posed on the basis of the existence of full first-order 
reductions not involving constraint mixing in the evolu- 
tion equations of the first-order variables. This intuitive 
line of reasoning, whether rigorously justified or not, nat- 
urally leads to the conclusion that the second-order ver- 
sion should also be ill-posed. Our results of Section [D] 
provide a rigorous basis for this intuitive inference. 

In the case of the BSSN equations the implications 
are less straightforward. In essence, the BSSN equations 
themselves constitute a partial first-order reduction of 
the ADM equations with constraint mixing in two senses. 
First, there is the use of a multiple of the Hamiltonian 
constraint in the evolution equation for the extrinsic cur- 
vature. Secondly, there is a use of the momentum con- 
straint in the evolution equations for the three new first- 
order variables T\ According to the latter play a 
critical role in the well-posed properties of the BSSN 
equations in a pseudospectral sense, whereas the for- 
mer is irrelevant to well-posedness in such a sense. This 
accounts for the difference in the well-posedness of the 
first-order reduction of BSSN as compared to ADM (but 
of course, the comparison is not completely fair because 
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BSSN is already a partial reduction). Yet, the constraint 
mixing in the evolution equations of the new first-order 
variables plays no role whatsoever in the well-posedness 
of the second-order problem for the three-metric com- 
ponents because it simply leads to a redundancy, as ex- 
plained in the main body of this article. Therefore, it 
is disappointing but hardly surprising that the second- 
order version of the BSSN equations does not turn out 
to be we ll p osed. In fact, this verifies the idea put for- 
ward in [lOj that the constraint mixing in the evolution 
of the first-order variables is crucial to the well-posedness 
of a first-order reduction of the Einstein equations. 

This idea is most strongly upheld by the results in the 
case of the generalized EC formulation and its second- 
order version. Here is a family of first-order reductions 
all of which are well posed for any multiple of the Hamil- 
tonian constraint added to the evolution equation for 
the extrinsic curvature. But they all have momentum- 
constraint mixing in the evolution equations for the 
new first-order variables fuj- Yet none of the second- 
order versions of these formulations turns out to be well 
posed. What the second-order versions are missing is 
precisely the evolution equations for the first-order vari- 
ables (which become redundant in the second-order ver- 



sion). For several reasons including the present ones, the 
constraints seem to be surfacing as key players in the 
initial value problem of the Einstein equations, contrary 
to their long-standing reputation as choosers of physical 
initial data but otherwise ignorable. 

Additionally, reflecting on the fact that the standard 
second-order version of the Einstein equations with har- 
monic slicing where all the components of the metric 
evolve according to wave equations is indeed well posed, 
one is led to conclude that the shift vector must play 
a crucial role in the well-posedness of second-order ver- 
sions of the 3+1 Einstein equations, in agreement with 
the intuition of the authors of [flJ where a dynamical 
shift choice is used to analyse some non-standard hyper- 
bolicity properties of the BSSN equations. This remains 
an open problem. 
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